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We derive an effective thermodynamic potential (ficff) at finite temperature (T 7^ 0) and zero 
quark-cliemical potential (pa = 0), using the singular-gauge instanton solution and Matsubara 
formula for A^c = 3 and Nj = 2 in the chiral limit. The momentum-dependent constituent-quark 
mass is also obtained as a function of T, employing the Harrington-Shepard caloron solution in 
the large- A'^c limit. In addition, we take into account the imaginary quark-chemical potential = 
A4, translated as the traced Polayakov-loop {^) as an order parameter for the Z{Nc) symmsetry, 
characterizing the confinement (intact) and deconfinement (spontaneously broken) phases. As a 
result, we observe the crossover of the chiral (x) order parameter and $. It also turns out that 
the critical temperature for the deconfinment phase transition, is lowered by about (5 — 10)% in 
comparison to the case with a constant constituent-quark mass. This behavior can be understood by 
considerable effects from the partial chiral restoration and nontrivial QCD vacuum on "I>. Numerical 
calculations show that the crossover transitions occur at (T^, T^) ~ (216, 227) MeV. 

PACS numbers: 12.38.Lg, 14.40.Aq 

Keywords: Thermodynamic potential, instanton, Polyakov loop 



I. INTRODUCTION 

The phase structure of QCD, as a function of temperature T and quark-chemical potential /u, represents the 
breaking patterns of the relevant symmetries in QCD. Simultaneously, each QCD phase can be characterized by the 
corresponding order parameters, reflecting the nature of the symmetries. In this sense, exploring the QCD phase 
diagram is of great importance in understanding strongly interacting systems. Especially, recent energetic progresses, 
achieved in the ultra-high energy experimental facilities, such as the RHIC, have triggered much interest to investigate 
the QCD phase structure in the vicinity of high T w Tc, whereas /i remains relatively small, being similar to the early 
universe. 

Starting from the first principle, the lattice QCD (LQCD) simulations must be a promising method to investigate 
this region (T ^ and ^ ~ 0) with less difficulties, such as the sign problem [THB]. Many attempts have been also 
done in various effective field-theoretical approaches [7Hl9]. Among them, interestingly enough, the Polyakov-loop- 
augmented Nambu-Jona-Lasinio (pNJL) model describes the crossover of the two different QCD order parameters 
for the chiral and Z(7Vc) symmetries, represented by the chiral condensate (qq) cx and the traced Polyakov loop 
{(/>) = $, respectively [Tl- fTTl [Hj . If the 'L{Nc) symmetry is intact, <I> becomes zero, indicating the confinement phase. 
On the contrary, provided that the symmetry is broken spontaneously, one has $ 7^ for the deconfinement one. 

Instanton model can be also thought as an appropriate framework to be employed for this finite-T subject, con- 
sidering that it has provided remarkable descriptions so far for various nonperturbative QCD and hadron properties. 
Note that the instanton solution at finite T, being periodic in Euclidean time, i.e., caloron turned out to be essential 
for this purpose [2Ql [21] . Nonetheless for its relevance, its practical application is still under development p2H27] . 
Confinement properties have been discussed as well with semi-classical objects, such as the meron (a half of regular- 
gauge instanton), by indicating the area law for the Wilson loop [551 [55] ■ The caloron with non-trivial holonomy, so 
called the Kraan-van Baal-Lee-Lu (KvBLL) caloron [20l[2T], was taken into account as a lump of dyons to understand 
the confinement [22 ] |23 ] [25] . 

In the present work, we want to construct an effective thermodynamic potential (r^off) finite T with /ip^ — 0, 
employing the instanton framework. Our strategy is rather simple and practical as follows: 

• Using the instanton distribution function at finite T from the caloron solution with trivial holonomy (the 
Harrington-Shepard caloron) [30[ 131] . we first compute the instanton density and average size of instanton 
as functions of T, resulting in that the instanton effect remains finite even beyond the critical temperature 
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Tc ^ Aqcd- Combining these ingredients, we finally obtain (A;(tliree momentum), T(temperature))-dependent 
constituent-quark mass M , M^^t^ which plays the most important role in the present approach. 

• In constructing ileff, we take into account a practical way, instead of using the caloron and its quark zero-mode 
solution: il^g is obtained by applying the Matsubara formula to the effective action, which is derived from the 
usual singular-gauge instanton solution at T = 0, as done usually in effective models [7 HTTI [Tll - fT^ . 

• The singular-gauge instanton solution is nothing to do with the confinement . On the contrary, it explains the 
nonperturbative QCD properties very well in terms of the spontaneous breakdown of chiral symmetry (SB^S). 
Hence, considering the chiral and 1{Nc) symmetries on the same footing as in the pNJL model, we introduce the 
imaginary quark-chemical potential /zi = A4, which corresponds to the uniform color gauge field in the Polyakov 
gauge. It will be translated later as the traced Polyakov loop as an order parameter for the spontaneous 
breakdown of Z(iVc) symmetry, i.e., the deconfinement phase transition. 

Considering all these ingredients, we can write a neat expression for D,cs as a function of T with M^^t- By solving 
the saddle-point equations with respect to the mean fields, we can draw the curves for and as functions of T. 
From the numerical calculations, we observe that the first-order deconfinment phase transition in pure-glue QCD is 
modified to the crossover one, according to the mixing of dynamical quarks and $ in ficff- In contrast, the mixing 
gives only negligible modification on a^. As a result, the crossover of the two different QCD order parameters is shown 
by the mixing. It also turns out that Tj, for the deconfinment phase transition (T^) is lowered by about 10% with 
Mk^T, in comparison to the case with a constant M (Afg.o ~ 350 MeV). From this observation, we can conclude that 
the nontrivial QCD vacuum contribution and partial chiral restoration play a considerable role even for deconfinment 
phase transition. The numerical results show that — 216 MeV and — 227 MeV. The discrepancy between them 
becomes about 10 MeV, which is rather larger than that computed in the local pNJL model [SHTU]. We note that the 
LQCD simulations provides smaller values than ours, whereas is comparable J]-^. We also find that cr^ depends 
much on the partial chiral restoration. 

We organize the present work as follows: We explain briefly the instanton properties and the effective chiral action 
derived at T = in Section II. The instanton distribution function at finite T is introduced, and M^^t is computed in 
Section III. In Section IV, utilizing Alk.T obtained in the previous Section, we construct an effective thermodynamic 
potential using the Matsubara formula with fj, = A4, which is translated into the traced Polyakov loop. The numerical 
results for and $ as functions of T and related discussions are given in Section V. The final Section is devoted to 
summary and conclusion. 



II. INSTANTON AT ZERO TEMPERATURE (T = 0) 

In this section, we want to explain the instanton framework for vacuum {T — and /iR — 0) bricfiy. First, we 
consider a grand canonical ensemble of dilute instantons, which indicate the tunneling in Euclidean time, as a dominant 
nonpertubative contribution of the QCD vacuum. In this framework, in principle, we have two phenomenolgical 
parameters such as the average instanton size po ~ 1/3 fni and inter-instanton distance ^0 ~ 1 [32| . where 
the subscript zero indicates vacuum. According to these values, the renomalization scale of this framework can be 
estimated as A sa 1 /po w 600 MeV. It is convenient to start with the quark zero mode in the presence of the instanton 
background in the chiral limit, assuming that the nonperturbative QCD properties are dominated by the zero mode: 

(z^-A//K°' = o, (1) 

where Ajj stands for a singular-gauge (anti)instanton solution, which can be obtained by minimizing the YM action, 
considering the self-dual equation and the SO (4) symmetry for the field strength tensor F^,^ in Euclidean space: 

where 77""^ indicates the 't Hooft symbol [35]. Provided that the zero mode saturates the low-energy hadron phenomena 

;(0) 



as mentioned, we can then write a single-quark propagator with the zero-mode solution, ^"1^? , approximately as 



1 *//^^// 

SiT = 1. — ~ So : , (3) 
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where m and 5*0 denote the small current-quark mass « and free-qaurk propagator Eq. ([s]) can be written 

alternatively, taking into account the Fourier transform of the zero-mode solution, as follows: 



S = 



1 



ip + iM{id) ' 



where the four momentum-dependent M (constituent-quark mass), M{id) is defined by 



(4) 



(5) 



where we have used the subscripts, representing three momentum and energy separately for convenience. Mq.o stands 
for M at zero momentum transfer at T = 0. It is worth mentioning that Mq q is a function of m in principle |33ji33]. 
Moreover, the meson-loop (ML) correction beyond the leading-iVc contributions turns out to be critical for the case 
with TO ^ |33l I35j . However, in the chiral limit as employed in the present work, Afo.o can be determined self- 
consistcntly without the ML correction, keeping only the leading-iVc contributions. The quark form factor, F{id) can 
be written in terms of the modified Bessel functions as follows: 



Fiid) = 2t 



Ioit)Ki{t) - h{t)Ko{t) 



t = 



Approximately, this form factor can be parametrized as follows: 

F{id) ^ 



(6) 



(7) 



Hereafter, we will make use of this paramterized form factor for numerical calculations. 

With all the ingredients taken into account so far, one can construct an effective partition function, which brings 
about the quark propagator given in Eq. (|3]) |32j : 



-2cff — I d\ D^p Dijj^ exp 



d'^x'4)\ip + im)ijj + X{Y+ + F") + iV(ln 



N 



XVM 



-1) 



(8) 



where we have introduced A as a Lagrangian multiplier to make the partition function exponent and an arbitrary 
massive parameter M to make the argument of the logarithm dimensionless. It is assumed that the numbers of the 
instanton and anti-instanton are the same, Nj ~ Nj = N/2 indicating the CP-invariant vacuum, so that N denotes 
the sum of the total number of the pseudo-particles. Since we are interested in the case for Nf = 2, the 2iV/-'t Hooft 
interaction can be casted into a four-quark instanton-induced interaction. From this Nambu-Jona-Lasinio (NJL) 
model-like four-quark interaction, one can construct an effective chiral action as a functional of quarks, and scalar 
and pseudo-scalar mesons, by an exact bosonization process |32| . After freezing out all the meson fields, except for 
the isosinglet scalar meson (ct), which is responsible for the spontaneous breakdown of chiral symmetry (SB^S), we 
have the following effective chiral action: 



1 - In 



N 
XVM 



2Vc7^ - 2VNrN 



J2nr 



(9) 



where SoS has been written in momentum space. More detailed explanations on the derivation of the effective action 
in Eq. Q from the partition function in Eq. dsl can be found in Refs. [311|331[3S]. Mq^q, as a function of the saddle- 
point value of a (ctq) and po, can be computea with the normalized instanton distribution function dN{p), which will 
be discussed in detail in the next section: 



Mo,o - Mao /" dpdN{p){2TTpf = Mcto(2^Po)' 



(10) 



Here, we have assumed that d^ip) ~ 5{p — po) in the large- iVc limit I36j . where M stands for a massive quantity, 
VX/ {2N^). Considering the saddle-point (self-consistent) equations with respect to A and a, one is led to the following 
relations: 



6X 
5a 







= 



N 
V 



d'^k 



^0,0 



F^{P 



cTo = 2A^c 



(27r)4 fc2 + Ml^F'^ik^) 
d^k Ml^F^jk^) 

(27r)4fc2 4.M2oP4(fc2)' 



no, 



(11) 



(12) 
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where we assigned the instanton number density (packing fraction) N/V as Uq. Considering the dilute-instanton 
ensemble, no « (200 MeV)'', po « 0.3 fm, and ~ 1-0 fm, and substituting Eq. ^ to Eq. ([ll|, we can determine 
Mg^o ~ 350 MeV self-consistently. Moreover, by equating Eqs. (Ill and (12), we have a relation 2ctq = uq at the 
saddle point. Using this relation, Eq. (10 1 can be rewritten as 



M, 



0.0 



/2^(27rpo)^ ~ 350 MeV. 



(13) 



Note that, in medium, ctq plays a role of an order parameter for the chiral phase transition as will be discussed in the 
following Sections, being similar to the chiral condensate (iq^q), which reads: 



(27r)4 + M^^F^ik^)' 



(14) 



III. INSTANTON AT FINITE TEMPERATURE (T / 0) 



Now, we want to discuss the instanton at finite T in the mean-field approach. In Ref. !12j, it was assumed that 
the relevant quantities such as no, po: and Mq.o remain the same with their vacuum values below and disappear 
beyond it. However, more realistic situation near Tc is quite complicated as indicated in Refs. |3 H 137 ) . in which the 
temporally periodic instantons with the trivial holonomy, i.e., the Harrington-Shepard caloron was employed [3U]. 
These works showed that the relevant quantities decrease with respect to T smoothly, and they beomce finite even 
beyond T,. 

We are in a position to consider the T-dependence of these relevant quantities, n, p, and M using the Harrington- 
Shepard caloron in detail. An instanton distribution function for arbitrary and Nf can be expressed with a 
Gaussian suppression factor as a function of T and an arbitrary instanton size p, employing the variational method 
for pure-glue QCD 

d{p, T) = Cn^ A^s P""' p'-' exp [~{AnT^ + hnp')p'] • (15) 
c 

Here, we considered the CP-invariant vacuum as before and assumed the same analytical form of the distribution 
function for both the insanton and anti-instanton. Note that n and p have been taken into account as functions of T 
implicitly. We also assigned the constant factor in the r.h.s. of the above equation as C for simplicity. The abbreviated 
notations are also given as: 

4 CQ „-1.68aRsAfc 

(3 = -51n[ARsPeut], P = -Hn[AR,s(i?)], Cn^_ = ^2<n-2)\(N - l)V ^^^^ 



An, = 



11 



^-Nr - 1 



27 

T 



b = 



llNr - 2N 



f 



(17) 



Note that we defined the one-loop inverse charge /3 and /3 at certain phenomenological cutoff values pcut and (i?) sa R. 
As will be shown, only /? is relevant in the following discussions and will be fixed self-consistently within the present 
framework. Ars stands for a scale, depending on a renormalization scheme, whereas V3 for the three-dimensional 
volume. Using the instanton distribution function in Eq. (15), we can compute the average value of the instanton 
size, p^ straightforwardly as follows [57] : 



Jdpp^d{p,T) 
Jdpd{p,T) 



AnT^ 



2/37n 



where v = [h — 4)/2. It can be easily shown that Eq. (18) satisfies the following asymptotic behaviors |37) : 



lim p^(T) 

T^O 



V 

f3jn 



hm p^T) 

T->oo 



AnT^' 



Note that the second relation of Eq. (19) provides a correct scale-temperature behavior at high T, i.e., l/p ■ 
Substituting Eq. (fTsl) into Eq. (IT5|, tf 



le distribution function can be evaluated further as: 



(18) 

(19) 
Acx T. 



dip, T)^C exp [- J-(T)p2] , F{T) = -A^T^ 



(20) 
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FIG. 1: n/no (solid) and p/po (dashed) as functions of AT for A^c = 3 (thick Unes) and 2 (thin lines). Here, we use no 
(0.2GeV)* and p= (0.6MeV)-\ 



The instanton density n can be computed self-consistently, using the following equation: 

n^F{T) = [CT{y)Y , 



(21) 



where we have replaced NT/V^ n, and r(t/) indicates the F-fucntion with an argument v. Note that C and f3 
can be determined easily using Eqs. (18 1 and (21 1, incorporating the vacuum values for n and p: C ~ 9.81 x 10~^ 
and /3 w 9.19. The numerical results for the normalized n/no (solid) and p/po (dashed) are shown as functions of 
T in Figure [l] for Nc = 2i (thick lines) and 2 (thin lines), separately. As shown in the figure, although n and p get 
decreased gradually as T increases, indicating the reduction of the instanton contributions, they remain finite even 
beyond ~ Aqcd « 200 MeV. 

Finally, in order to estimate the T-dependence of M, it is necessary to consider the normalized distribution function 
(see Eq. (lOl), defined as follows: 



dN{p,T) = 



d{p,T) 
Jdpd{p,T) 



P 



h-5 



J-'^(T)exp[-J-(T)p2] 



r(^) 



(22) 



The numerical result for djq{p^ T) is given in Figure [2j As shown in the figure, the peak of the dj^{p, T) surface moves 
to smaller p as T increases on the p-T plane. Now, we want to employ the large- A^c limit to simplify the expression of 
d^iPi T) as done for the T = case in the previous Section. Since the parameter b is in the order of 0{Nc) as shown 
in Eq. (161, it becomes infinity as Nc — >■ oo, and the same for i/. In this limit, as understood from Eq. (|22[), d^iPjT) 



can be approximated as a (5-function 



lim 



dN{p,T) = 6{p--p{T)). 



(23) 



The trajectory of this ^-function projected on the p-T plane is depicetd in Figure[2]by the thick-solid line. Substituting 
Eq. (23) into Eq. (10) and using Eq. ([7|, we can write the (fc, T)-dependent M as follows: 



Mi 



h.kA 



Mo.o 



'n p 



2 + (fc2 + kl) p2 



(24) 



Mo- 



where we use Mq q ~ 350 MeV, considering that Eq. (24) must be identical to Eq. (jsl as T — > 0. Moreover, we 
assumed that the Lagrangian multiplier A, equivalently M, is independent on T. Note that the fourth component of 
the momentum, ^4 will be replaced as a function of T by employing the Matsubara frequency in the next Section. 

Here, we discuss a caveat in our model. In deriving the momentum- and temperature-dependent constituent-quark 
mass in Eq. (24), we have used a fact that the Harrington- Shepard caloron is the zero-mode solution of the Dirac 
equation with the trivial holonomy in the presence of the instanton background at finite T. However, as mentioned 
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FIG. 2: dN{p,T) on the p ~T plane. The trajectory on the curve represents p{T) in Eq. (181 



in Section I, we will employ the Polyakov loop in Section IV, corresponding to the nontrivial holonomy, resulting in 
that the Harrington-Shepard caloron is not the exact zero-mode solution of the Dirac equation anymore. Instead, the 
KvBLL caloron becomes the exact one. 

Since the remedy for this problematic issue is quite complicated and beyond our scope in the present work, we 
assume that the Harrington-Shepard caloron is almost the zero-mode solution even with the nontrivial holonomy so 
that Eq. (24) is qualitatively apphcable. 



IV. EFFECTIVE THERMODYNAMIC POTENTIAL 

In this section, we construct an effective thermodynamic potential f2eff, considering all the ingredients discussed so 
far. In addition, we take into account the imaginary quark chemical potential (/ii = A4), which will be translated 
as the traced Polyakov loop as an order parameter for the deconfinment phase transition [71 [T3]- Moreover, this 
corresponds to the uniform color gauge field, induced in the Polyakov gauge, = ((3,^4). All calculations will be 
performed in the case for Nc = 3, Nf = 2, and /iR = in the chiral limit for the leading-iVc contributions. In order to 
evaluate f2cff as a function of T from the effective action in Eq. ([9|, we employ the anti-periodic Matsubara formula 
for fermions in Euclidean space: 



(2^ 



(25) 



where the Matsubara frequency reads — (2n -I- \)i:T . We also include an imaginary quark-chemical potential |71 
[TU] . which can be identified as the fourth component of SU(A^c) gauge field (A4) in Euclidean space, resulting in 
a replacement fc — > fc — ^4 in rieff- Applying these ingredients to Eq. (|9]), we can have the following effective 
thermodynamic potential per unit volume with the mixing of the dynamical quarks and A^: 



q-|-A4 
eff 



2o-^ - iV/T 



E 



(2^ 



Tr^ In 



(fc - A^f -f Ml; 

2^2 



(26) 



where cr is a function of T as well. Note that we have ignored Ai^ in M for simplicity. After a straightforward 
manipulation, using the relations and identities given in Appendix, one is led to the following expression: 



2a' 



d^k 



rTr, 



k.T 



T 



In 



1 



1 



(27) 



Here, we have used that 



k,T 



1/2 



M, 



k,T 



M, 



0,T 



(28) 



2 + fc2p2_ 

e nun 

k\. From the figure, it turns out that the k dependence of M becomes weak as T 



where Mq t is defined in Eq. (24|). For more details, see also Appendix. The numerical result for M^ t is shown in 
Figure |3] as a function of T and 



FIG. 3: Mk,T in Eq. 1 28 1 as a function of T and the absolute alue of three momentum |fc| [GeV]. 



increases. At T — 0.4 GeV, the k dependence does not appear at all. This behavior can be understood as follows: 
The k dependence of M is generated from the quark-instanton interaction, i.e, the delocalization of the quark fields in 
the instanton ensemble. As the instanton ensemble becomes more dilute as T increases (smaller p, equivalently) , the 
interaction probability decreases, resulting in the reduction of the k dependence. However, note that the dependence 
still remains visible around T^, ~ Aqcd ss 200 MeV as shown in the figure. 

Now, we are in a position to consider the traced Polyakov loop 0, defined in a SU(7Vc) gauge group as: 



-Tr.exp — 



-Tr.exp — 



(29) 



Taking into account the Polyakov gauge, A4 is diagonal in a A^c x Nc matrix. For instance, the perturbative YM 
potential can be expressed by this physical quantity [31] and prefers the deconfinment phase [33110]: (</>) becomes 
finite according to the spontaneous breakdown of the Zn^ symmetry with the trivial holonomy. In contrast, if the 
symmetry is intact, one finds ((/)) = 0, indicating the confining phase with the non-trivial holonomy. 

Thus, (j) plays the role of an exact order parameter for the 1,{Nc) symmetry for pure-glue QCD, in which the quark 
degree of freedom is decoupled according to its infinitely heavy mass. It is worth mentioning that the 1,(Nc) symmetry 
is broken explicitly in the presence of dynamical quarks with finite mass, considering the anti-symmetric nature of 
fcrmions, resulting in that (0) is not an exact order parameter for Zj^^ symmetry any more. However, from the 
phenomenological point of view, incorporating (j) and dynamical quarks has been quite successful to a certain extent 
to explain various features of the QCD phases transition: The crossover near Tc for instance for Nf = 2. Hence, in 
the present work as done in the pNJL model, we have incorporated the instanton-based model with the SB^S and 
the traced Polyakov loop as an order parameter for the 'Z(Nc) symmetry. The trace over color space in Eq. (27) can 
be evaluated further in terms of and cj)* using Eq. ( 29 ) as follows: 



Tr.ln 



In 



-In 



I + A^e U* + 



1 + Nc [ (j) + (t)* e- 



(30) 



On top of Jlgf^^" (now becoming Q'^^'^), an additional pure-glue effective thermodynamic potential was suggested 
in Refs. [gHTTl as a function of and 0*: 



_ _rji4: 

"cff — ^ 



&2(r), 



(31) 



where the coefficient 69 is a function of T: 



62 (T) = ao +ai 



'To' 




'% 


1 




To" 


T 


+ a2 


T 




V as 


T 



(32) 
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Here, Tq denotes critical T in pure-glue QCD, resulting in Tq = 270 MeV at which the first-order phase transition 
occurs. Note that Tq is different from Tc, which will be computed later in the presence of the mixing of the dynamical 
quarks and (j). The coefficients, a and b, are listed in Table [l] [H]. This parameterization of the effective potential in 
Eq. (31) bears the Zn^ symmetry, conserved in pure-glue QCD by construction, and works qualitatively well up to 



T « (2 — 3) Tc, from which the transverse gluons come into play considerably. 

Finally, substituting Eq. (30) into Eq. ( p7| ) and adding Eq. (31) to it, we arrive at the following expression for 
the effective thermodynamic potential with the two order parameters, for the chiral phase and (j> (0*) for the 
deconfinment phase transitions, at finite T and /ir = 0: 



Nr. 



(Pk 

(2^ 



k.T 



T 



T 



(27r)3 

d^k 
(2^ 

b2{T) 



In 



In 



1 + iVc * + * 



($$) 



iVc + * 
6 ^ 



$3 



4 ^ ' 



(33) 



where we have replaced </> into its mean value {(j>) = $ as done in Rcfs. [SUlOj. Although, this expression for ilcs is 
very similar to those given in Refs. [8 -ITIIII], ours is distinctive from them quantitatively in several points; 

• The scale parameter of the model A « 1/p is obtained as a function of T by solving the instanton distribution 
function as discussed in the previous Section. Moreover, M is expressed as a function of k and T {M^^t), rather 
than a constant, manifesting the partial chiral restoration and nontrivial QCD vacuum contributions. 



• Consequently, there appears no divergence in the energy integral cx J cPk E^^t in Eq. (33) by virtue of the 
fc-dependent M, which plays the role of an intrinsic ultraviolet (UV) regulator, being different from other 
local-interaction models, such as the usual pNJL model. 

• At the smae time, as for the 2Nf-'t Hooft interaction, the quark-meson coupling strength depends on k as well 
as T, being different from that in other models, in which it is a constant value fixed at zero T. 

• All the relevant quantities at zero T, Mq^q, ctq and uq for instance, are determined self-consistently by solving 
the saddle-point equations, Eqs. (11) and (12). 



Now, we evalute the equations of motion with respect to the mean fields, a, $, and $, by minimizing flcS' 



0, 



0, 



0. 



(34) 



(5$ ' ,5$ 5c7 

From the first two equations, it can be easily seen that $ = $ at the saddle point, and we have the relation as follows: 

4:NaNf (e- 



[64 $^ - &3 $2 - b2{T) = y 



d^k 



{2n) 



1 + Nc<^> 1 



(35) 



Similarly, the last equation in Eq. (34) gives 



2cr^ = N^Nf 



d^k Ml^ 
(27r)3 Ek^T 



2 $e- 



1 + Nc^ 1 



2$e- 



I 1 + e T I e T + I 



(36) 





ai 


02 


a4 


63 


64 


6.75 


-1.95 


2.63 


-7.44 


0.75 


7.50 



TABLE I: Coefficients for f^f^ in Eqs. (31l and (32l, taken from Refs. [8HT0). 
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FIG. 4: <I> (increasing curves) and normalizaed (decreasing ones) for To = 270 MeV using Mh.T- The solid and dashed lines 
indicate the cases with and without the dynamical quark and $ mixing, respectively. 



where we have used the following relation 



da 



(37) 



considering Eqs. (10) and (13). From Eq. (36), one can immediately see the following asymptotic behaviors: 



lim [2cr2l 

T-fO 



(2^)3 Ek,T 



as $ 0, 



lim \2a^] = as $ -> 1, 

T-i-oo 



(38) 



showing appropriate chiral properties, as expected. 



V. NUMERICAL RESULTS AND DISCUSSIONS 



In this Section, we present the numerical results for the two different order parameters, and $ as functions of T, 
and related discussions are given as well. Hereafter, we will used the normalized value for ct^, cr^/cg for convenience. 
Firstly, in Figure |4j we draw them for the cases with (solid lines) and without (dashed lines) the mixing of the 
dynamical quark and $, using Mfc.T- Here, we employ To — 270 MeV for the pure-glue potential in Eq. (31). In 



order for the unmixed case (dashed lines) , in whic h th e quarks are decoupled from the pure gluodynamics, we set M 
infinite, equivalently E^.t — > oo in Eqs. (35) and (36). As shown in the figure, even with or without the mixing, cr^ 



showed the crossover for the chiral phase transition. On the contrary, the phase transition pattern for $ turns out 
to be different, depending on the mixing. In this sense, the crossover of $ is caused by the mixing, as suggested by 
the pN JL modelsTTHTUl and shown by the LQCD simulations with dynamical quarks [5] . An interesting behavior 
shown in Figure |4 is that the mixing effect increases then decreases for as indicated in the solid and dashed lines. 
This tendency comes from the combination of the opposite behaviors of 4> (increasing) and e~^'='^/"^ (decreasing) in 
cr2 in Eq. ([set. 

In Figure|5 we present the numerical results for $ as a function of T for different types of M, as listed in Table [11} We 
again employed Tq = 270 MeV. $ for the pure-glue potential, showing the first-order deconfinment phase transition. 



is also given in Figure [5] for comparison. From the figure, it is clearly shown that the T- and/or fc-dependent M make 
$ shifted to lower T, resulting in lowering T^ (we will discuss the numerical values for Tc later in detail). 

In order to see the effects of {k,T) dependence in M, we take a careful look on Eq. (35). For simplicity, we 
approximate it considering only the leading contributions as follows: 



T^ [64 - 63 $^ - 62 (r) $] sa m^Nf 



(fik 
(2^ 



-Eu.t/T 



(39) 
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As known from Figure [s] and Eq. (28), when we take into account the T- and/or fc-dependent M, the strength of M 
decreases as fe and/or T increases, and the same for Ekj-. As a result, the integrand in the r.h.s. of Eq. (39) tends 
to be larger, as k and/or T increases, than that with a constant M, Mq q « 350 MeV. To make things clear, we put 
T = To in Eq. ( 39 1 for example, then have 



/(Pk f 



(27r)3 



--E'o,o/Tb 



(40) 



where we have used the notation, E^f^ — Mqq + k^. In order to satisfy the relations of Eq. (40), $ for A, must 
be bigger than aX T — Tq, since <i> is positive definite and 64 ^ 63. This observation is also true for arbitrary T, 
except for the limiting cases, T — ^ or 00. Ats the same time, therefore, this situation can be understood as that $ 
is shifted downward almost horizontally: Tc is lowered consequently. 

However, this lowering Tc behavior must be understood separately for the k and T dependences in M, since the 
both curves with Mq^t and Mk^ show it, as depicted in Figure [Sj 

• T dependence: Afo,T 

It tells us that the partial chiral restoration, indicated by decreasing M oc with respect to T, has effects on 
the deconfinment phase transition, although these two phase-transition mechanisms are believed to be different 
to each other. Here is a microscopic explanation for this: As T increases, the pseudo-particle (instanton) become 
smaller in its size ~ p, decreasing QCD vacuum contribution simultaneously, resulting in that the quarks are 
enable to travel more freely with less interactions with the instantons, and strings attached to each quark are 
extended more at lower string tension, in comparison to the case without the partial chiral restoration in M. 
Hence, the condensation of the strings can happen easily at lower T, toward the deconfinment phase. 

• k dependence: Mfe^o 

The k dependence in M is originated from the delocalization of quark fields in the presence of the instanton 
background [32], not from T-related effects. As momentum transfer increases, quarks become lighter c>c l/fe"*, 
loosing the nontrivial QCD vacuum contributions. In other words, the instanton effect is weakened seemingly in 
the k integrals. Consequently, being similar to the T-dependence case, the deconfinment phase transition takes 
place at lower T. 

From these explanations, one is led to a conclusion that the QCD vacuum contributions and the partial chiral 
restoration play considerable roles for the deconfinment phase transition to a certain extent in the presence of the 
mixing. 

Now, we are in a position to calculate the critical temperatures for the crossover transitions numerically. They 
are determined by the chiral and Polyakov susceptibilities as in Ref. j7J. Being almost equivalently, they can be also 
obtained from the maximum values of d^/dT and da^ /dT for the deconfine ment and chiral phases, respectively as in 

we list them for <I> and cr^, assigned 



III 



Refs. [8l-[10l 144] ■ In the present work, we employ the later method. In Table 
as and T^, respectively, for each type of M. Note that we do not show the numerical results for for the cases 
with Mo,T and Afo,o, since they are UV divergent, proportional to / k^dk, unless a cutoff is introduced by hand. 

As discussed previously, from the table, one can see clearly that Tc is lowered by inclusion of the (fe, T) dependence 
in M. It turns out that the shift of T^ is about 10%, (240 227) MeV for Mo,o Mj. t- Interestingly, if we take into 
account full (fe, T) dependence, and T^ get closer to each other as shown in Table |lll| (T*, T'"') = (227, 216) MeV. 
About 5% discrepancy (^ 10 MeV) is observed between them, showing a tendency T^ < Tf. By turning off the T 
dependence in M {A'Ik,o), the discrepancy between T^ and becomes larger up to about 15%, {T'^,T'^ ) = (225, 265) 
MeV. For interpreting this behavior, we take a look on the T dependence of for the cases with M^^t and M^^, and 
draw the numerical results in Figure |6] For comparison, we also draw $ for the two cases. Note that cr^ show obvious 
difference for the cases with M^^t and Mk^. This is the reason why the discrepancy between and T^ is so large 
for Mkfi. At the same time, this strong dependence on T for the chiral phase transition interprets the larger shift of 
T^, (265 — ?► 216) MeV. From this observation, we can conclude that the T dependence in M plays an important role 
in exploring chiral phase transition at finite T. 



Mk,o 



Mo- 



Mo 



Mo,o \^^] 



350 MeV 



TABLE II: Notations for M given in Figure [5] 
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FIG. 5: $ as a function of T for To = 270 MeV, employing Affej (solid), Mk^ (dotted), Mo,t (dashed), and Mo,o (long-dashed). 
We also draw "I> for the pure-glue case (dot-dashed). 




FIG. 6: $ (increasing curves) and (J'^/cto (decreasing ones) as functions of T for To = 270 MeV, employing Mh.T (solid) and 
Affc,o (long-dashed). 



It is worth mentioning that, from the LQCD analyses, it turned out that ~ 180 MeV for Nj — 2 using the 
clover-improved Willson fermions [3]. Also, using the renormalization-group (RG) improved action, it was found that 
« 171 MeV fS]. These values are significantly smaller than ours by (10 — 20)%. If this is the case, one may need 
more strong T dependence for M in the present approach, since, taking a look on Eq. (36), the T dependence of 



tj is mainly governed by the behavior of M. In other words, the instanton effects must decrease much faster as T 
increases. There can be several possible scenarios to satisfy this condition: 



• If we consider a correct instanton distribution function, rather than the simplified one in Eq 
the large- A^c hmit, may be lowered. To test this, we have the ratio of M computed with t 



( 23 ) according to 
le correct one in 



To = 270 MeV 






Mo,T 


Mo,o 


$ 


= 227 MeV 


= 225 MeV 


T<? = 230 MeV 


Tf = 240 MeV 


a /(Jo 


rf = 216 MeV 


rf = 265 MeV 







TABLE III: Tc computed from <1> and a^/a'^ for To = 270 MeV. 
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Eq. (22) and Mo,t in Eq. (24): 



M, 



0,T 



dp p e 



(0.46- 0.47). 



(41) 



We, however, verify that this modification does not work for lowering T^^ whereas the strength of M is reduced 
approximately by half. 

• Additional T-dependent terms can be taken into account. For instance, we assumed that the Lgrangian multiplier 
A in Eq. ([8| is independent on T in deriving Affc,T in Eq. (24). If it has the T dependence, we can modify Affc,T 
as follows: 



(42) 



Quantitatively, the T-dependent A can not be determined self-consistently within the framework. From a very 
rough estimation, based on the assumption that hq in Eq. (11) converges with a brute cutoff Aw 1/poi we can 
obtain a relation Aq oc 1/pq |32j . Using this assumption, the expression for M^^t in Eq. (24) is modified into 



'0 



/n p 



Mo.o 



/np 



no po 



(43) 



It turns out that this rough assumption makes things worse: is shifted to a larger value as expected in 
Eq. (43). In addition, it is also possible that the reduced coefficient C in Eq. (15) and the one-loop inverse 



charge (3 in Eq. ( 16 1 are functions of T, although it is beyond our scope in the present work. 



We may not ignore the Matsubara frequency Wn in the denominator in the l.h.s. of Eq. (44 1 in Appendix 



In the presence of the mixing, $ may provide effects on the instanton distribution function, becoming an 
exponentially decreasing function, not a gaussian one. We, however, do not perform quantitative calculations 
for these possibilities here and leave them for future works. 



• Additional T dependence can be added to the instanton distribution function in Eq. ( 15 ) according to the 
fcrmion overlap matrix, if one considers full QCD in computing the distribution function [42] . being different 
from the present work based on the variational method in pure-glue QCD (the Harrington-Shepard caloron) |31j . 
Moreover, the instanton clustering, which was suggested as a main contribution for the chiral phase transition |42l 
I43j and not taken into account here, may be responsible for lowering T^. 

• We note that Tq for the pure-glue potential can be chosen as a smaller value than 270 MeV, which has been used 
throughout in the present work. As shown in Refs. [5HTU]. by taking Tq = 190 MeV, the computed values for 
became compatible with those from the LQCD simulations. Although we have not presented detailed results 
for lower Tg, we could obtain = 194 MeV using M^.t at Tq — 200 MeV, showing about 10% decreasing. 

As for estimated in the LQCD, using the clover-improved Wilson fermions similarly, it was determined about 210 
MeV [5], which is rather compatible with ours. It was observed in the usual local pNJL model [51-[TU] that Tc taken 
from the two different order parameters are almost consistent: K,Tf k, 220 MeV for Tq — 270 MeV. It turned out 
that Tf = 215 MeV in Ref. [IJ with pNJL. 

Finally, we compare our results for $ with the LQCD data from Refs. j45j and [15] in Figure [7j in which a full and 
quenched calculations were done for Nf ~ 2. In their works, it was observed that = 202 MeV, which is about 10% 
lower than ours, 227 MeV for Nf = 2, and 270 MeV for Nf = 0. As shown in the figure, the LQCD data, indicated 
by □ (full) and A (quenched), are in qualitative agreement with the present results, but not quantitative. Especially, 
our result for the mixed case deviates much from it for the region T > Tq. 



VI. SUMMARY AND CONCLUSION 



In the present work, we have attempted to construct an effective thermodynamic potential il^s at finite T and 
zero quark-chemical potential (/iR = 0) in the chiral limit. We restricted ourselves to A^c = 3 and Nf = 2. Motivated 
by the Polyakov-loop-augmented Nambu-Jona-Lasinio (pNJL) model, we wanted to incorporate two different order 
parameters, cr^ and $, which characterize the chiral and deconfinment phase transitions, respectively. 

In order to discuss the spontaneous breakdown of chiral symmetry at finite T, we employed the singular-gauge 
instanton solution, and the fermionic Matsubara formula to express the effective chiral action as a function of T. We 
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1.0 



11 



□q 



0.75 



1.00 



1.25 



1.50 



FIG. 7: "1> as a function of T for To = 270. The solid and dashed Hnes indicate the cases with and without the dynamical 
quark and $ mixing, respectively. The notations □ and A indicate the full and quenched lattice data, respectively, taken from 
Refs. [45] and [46]. 

employed the instanton-distribution function, derived from the Harrington-Shepard caloron, to obtain the instanton 
density and average instanton size as functions of T. It turned out that these two quantities decreased but finite, 
indicating that the instanton effect survives even beyond |3T|. The (fe, T)-dependent M, Mk^T was derived in the 
large- iVc limit. We found that the k dependence of M becomes weaker as T increases. At the same time, the absolute 
value of Mk,T was also reduced with respect to T. To include the Polyakov loop as an order parameter for the Z{Nc) 
symmetry, we took into account imaginary quark-chemical potential fij = A^, which was translated as the traced 
Polyakov loop Combining all these ingredients, we could construct Jlcff with an additional pure-glue SU(iVc) gauge 
potential. By minimizing fieS with respect to external fields such as a and we could compute cr^ and $ as functions 
of T numerically. From the various numerical results we have found the followings: 

• In the presence of the mixing of the dynamical quarks and $, we observed that <& is very sensitive to the mixing, 
showing the crosssover and first-order transitions with and without it, respectively. In contrast, is insensitive 
to it, indicating the crossover phase transition. 

• M was expressed as a decreasing function of T as well as k. Due to this, was lowered by about (5 — 10)%, 
in comparison to that with constant mass Afo,o ~ 350 MeV. From this observation, we explain this lowering 

by that the nontrivial QCD vacuum contributions and partial chiral restoration play a significant role even in 
the deconfinment phase transition. 

• If the (fe, T) dependence had been fully taken into account, we found that = 227 MeV and = 216 MeV. 
The discrepancy between them became about 10 MeV, which was rather larger than that computed in the pNJL 
model. We also note that the LQCD simulations presented smaller than ours, whereas was compatible. 

• Finally, we observed that was depending much on T. Again, the partial chiral restoration turned out to be 
crucial to make proper results for the chiral phase transition. 

Consequently, from the present work, we could learn that the partial chiral restoration and nontrivial QCD vac- 
uum effects must be taken into account appropriately to investigate the breaking patterns of the chiral and Z{Nc) 
symmetries at finite T. As a next step, we attempt to include the finite quark-chemical potential (/l(r), giving a full 
description for the QCD phase diagram on the /x^-T plane. In addition to it, the finite current-quark mass, m is 
also under consideration beyond the chiral limit. This is important, since the explicit breakdown of the flavor SU(3) 
symmetry modifies the QCD phase diagram to a good extent. However, including finite m into the present framework 
has a huddle: One needs to consider the meson-loop corrections, which make significant modification on the physical 
quantities such as the chiral susceptibility, in comparison to those in the chiral limit [3 5) . Related works are underway 
and appear elsewhere. 
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Appendix 

The logarithm in Eq. (26) can be evaluated using the parametrized quark form factor, given in Eq. Q as follows: 



In 



2 + (A;2 + iD2)p2 



In 



2 + fe2p2 



(44) 



We verified that difference between the r.h.s. and l.h.s. of Eq. (44) is almost negligible. We used the following 
notations as 



0,0 



/n p 



(45) 



where /i stands for complex quark-chemical potential, /i <C 1. In the present case, we choose /x — lA^. In deriving 
Eq. (26) from Eq. (27 1, we also used the following Matsubara-sum indentity: 



T ^ In 



2 

k,T 



rp2 



Ek.T + T 



ln[l + e T ] + ln[l + ^ ] 



where the energy of a quark is given by 



0,T 



2 + fc2p2 



1/2 



(46) 
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